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ABSTRACT 


This  study  examines  an  alternative  method  for  modeling  in-tube  condensation  heat  transfer.  The 
method’s  fundamental  premise  is  that  the  heat  flux  is  directly  proportional  to  the  difference 
between  the  enthalpy  of  the  vapor  at  saturated  free  stream  conditions  and  the  enthalpy  of  the 
condensed  subcooled  liquid  at  the  wall.  The  proportionality  coefficient  is  a pseudo  condensation 
mass  velocity.  The  traditional  method  for  modeling  in-tube  condensation  uses  Newton’s  law  of 
cooling,  which  defines  the  heat  flux  to  be  directly  proportional  to  the  temperature  difference 
rather  than  the  enthalpy  difference.  Kraay  and  Rite  have  used  an  enthalpy-based  expression  for 
the  condensation  heat  flux  where  the  condensation  mass  velocity  was  derived  from  kinetic 
theory.  The  present  research  demonstrates  that  the  kinetic  theory  approach  does  not  yield 
realistic  results.  Consequently,  this  study  presents  an  alternate  energy  balance  approach  is 
subsequently  used  to  solve  for  an  enthalpy  based  condensation  heat  flux.  The  method  reveals 
that  the  heat  transfer  process  depends  strongly  on  the  derivative  of  the  vapor  quality  with  respect 
to  the  distance  along  the  tube  axis.  Based  on  this  derivation,  a correlation  is  presented  that  is 
simpler  than  those  derived  from  traditional  methods.  The  temperature  difference  between  the 
free  stream  vapor  and  the  liquid  at  the  wall  does  not  appear  in  the  correlation. 

Keywords:  in-tube  condensation,  kinetic  theory,  enthalpy-based  method,  correlation. 
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NOMENCLATURE 


English  symbols 


A 

cp 

Dh 

G 

i 

*fg 

h 

h,2(p 

Ja 

Jm 

k 

L 

m 

M 

Nf 

Nu 

P 

Pe 

Pr 


pw 


area  (m2) 

specific  heat  (J  kg'1  K'1) 
hydraulic  diameter  (m) 
mass  velocity  (kg  m'2  s'1) 
enthalpy  (J  kg’1) 

latent  heat  of  condensation  (J  kg'1) 

local  single-phase  heat-transfer  coefficient  (W  m'2  K'1) 

local  two-phase  heat-transfer  coefficient  (W  m'2  K1) 

if 

refrigerant  Jakob  number (-) 

cp(Tsat-TJ 

mass  flux  per  unit  time  and  per  unit  area  (kg  m'2  s'1) 

thermal  conductivity  (W  m'1  K'1) 

tube  length  (m) 

mass  flow  rate  (kg  s'1) 

molecular  weight  (kg  kmol'1) 

number  of  fins  (-) 

local  Nusselt  number  based  on  Dh 

pressure  (Pa)  (-) 

Peclet  number  (Re-  Pr)(-) 


liquid  refrigerant  Prandtl  number 


cp.iPi 
k 1 


wetted  perimeter  (m) 
local  heat  flux  (W  m'  ) 
inside  radius  of  the  pipe  (m) 

radius  of  the  cylindrical  surface  equivalent  to  the  actual  liquid-vapor  interface  (m) 


R 

Re 

S 

Sv 

St 

T 

u°° 

V 

x 

z 


universal  gas  constant  (N  m mole'1  K1) 

GDh 

liquid  refrigerant  Reynolds  number (-) 

Pi 

perimeter  of  one  fin  and  channel  - measured  perpendicular  to  the  fin  axis  (m) 

v,  - vf 

non-dimensional  refrigerant  specific  volume (-) 

h 

Stanton  number, (-) 

PVcP 

temperature  (K) 

velocity  of  the  vapor  at  the  pipe  entrance  (ms'1) 
generic  velocity  (ms1) 
thermodynamic  mass  quality  (-) 
axial  distance  (m) 
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Greek  symbols 

T coefficient  introduced  in  Eq.  A.4  (-) 
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A 

X 

V 

V 

P 

a 

turbulent  liquid  thickness  (m) 

X-8  (m) 

thickness  of  the  liquid  layer  (m) 
dynamic  viscosity  (kg  m'1  s'1) 

equilibrium  two-phase  specific  volume  xvg+(l-x)vi  (m  kg'  ) 
density  (kg  m'3) 

condensation,  or  evaporation  coefficient  introduced  in  Eq.  A.4  (-) 

Subscripts 

c 

e 

f 

a 

& 

glide 

ig 

1 

lat 

r 

ref 

sat 

sens 

tot 

w 

condensation 

evaporation 

saturated  liquid  phase 

saturated  vapor  phase 

due  to  temperature  glide 

ideal  gas 

subcooled  liquid 

latent 

refrigerant 

reference 

saturated 

sensible 

total  (liquid  and  vapor) 
wall 
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Introduction 


The  traditional  method  to  modeling  two-phase  heat  flux  was  proposed  in  the  early  1700s  in  an 
anonymous  paper  generally  attributed  to  Isaac  Newton  (Newton,  1701).  Newton  describes  an 
experiment  where  the  heat  transfer  between  two  fluids  at  different  temperatures  is  observed  to  be 
proportional  to  their  temperature  difference.  Newton’s  law  of  cooling  is  the  mathematical 
expression  of  this  observation.  The  Newton’s  law  of  cooling  for  a condensation  process  dictates 
that  the  heat  flux  released  by  the  condensing  fluid  ( q" ) can  be  expressed  as  the  product  of  the 
temperature  difference  between  the  saturated  fluid  and  the  wall  ( TsarTw ) and  a coefficient  of 
proportionality,  which  is  called  heat  transfer  coefficient,  i.e.: 

q"=K(Tsa:~TJ  (1) 

The  traditional  approach  used  in  studying  the  two-phase  convective  heat  transfer  process  consists 
of  predicting  the  heat  transfer  coefficient  (h2(p),  by  correlating  it  to  dimensionless  parameters.  The 
two-phase  heat  transfer  coefficient  is  a function  of  several  parameters,  including  vapor  quality, 
mass  flow-rate,  pressure,  temperature,  and  the  thermodynamic  and  transport  properties  of  the 
fluid.  Additionally,  the  form  of  the  correlation  varies  with  the  two-phase  flow  regimes  and  the 
theoretical,  semi-theoretical  or  empirical  approach  used  in  modeling  the  heat  transfer.  According 
to  Hinde  et  al.  (1992),  different  forms  of  correlations  can  produce  predictions  that  differ  by  a 
factor  of  two  or  three.  A comprehensive  review  of  the  state  of  the  art  is  presented  in  Dobson  and 
Chato  (1998).  The  drawback  to  the  traditional  approach  is  that  the  heat  transfer  coefficient 
correlation  is  affected  by  (1)  the  choice  of  the  dimensionless  parameters  for  its  functional  form 
and  (2)  the  choice  of  the  approach  to  be  used.  It  is  clear  that  the  considerable  difference  between 
prediction  methods  causes  uncertainty  in  the  rating  and  sizing  of  heat  exchangers.  As  a result, 
condensers  are  usually  oversized  to  meet  the  required  duty. 

To  overcome  the  disadvantages  of  the  traditional  approach  to  correlating  the  two-phase  heat 
transfer  process,  Kraay  and  Rite  (1997)  proposed  a new  methodology  in  which  the  heat  flux  is 
assumed  to  be  proportional  to  the  difference  between  the  enthalpy  of  the  vapor  at  free  stream 
conditions  and  the  enthalpy  of  the  condensed  liquid  at  the  wall: 

(2) 
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The  coefficient  of  proportionality,  Jmi  represents  the  flux  of  mass  per  unit  time  or  mass  velocity 
from  a point  at  the  free  stream  conditions  to  a point  at  wall  conditions.  The  advantage  of  the 
Kraay  and  Ray  method  for  modeling  condensation  is  that  it  accounts  separately  for  the  phase 
change  heat  transfer  and  the  liquid  subcooling  heat  transfer.  The  traditional  temperature  based 
heat  transfer  equation  cannot  be  used  to  describe  the  phase  change  alone  because  the  phase 
change  occurs  at  essentially  no  temperature  difference. 

Kraay  and  Rite  (1997)  claim  that  their  approach  allows  the  parameter  that  they  refer  to  as  a “new 
heat  transfer  coefficient,”  Jm,  to  depend  only  on  flow  rate,  fluid  properties,  and  wall  geometry, 
rather  that  being  a function  of  local  vapor  quality.  They  studied  the  latent  heat  transfer  at  the 
interface  by  using  the  kinetic  theory  of  condensation  (this  approach  is  described  in  the 
Appendices).  It  is  difficult  to  rigorously  apply  the  kinetic  theory  to  convective  condensation.  For 
example,  the  kinetic  theory  of  condensation  is  well  developed  for  the  cases  in  which  the  velocity 
of  the  condensing  vapor  has  no  component  tangent  to  the  liquid-vapor  interface  and  where  the 
condensation  process  is  independent  of  vapor  quality.  Unfortunately,  in-tube  condensation  is  a 
strong  function  of  vapor  quality  and  exhibits  significant  tangential  vapor  velocities.  To 
compensate  for  this,  Kraay  and  Rite  (1997)  revert  to  several  assumptions  to  adapt  the  kinetic 
theory  to  in-tube  condensation,  which  are  not  always  easy  to  justify. 

The  current  study  has  two  objectives.  The  first  objective  is  to  review  the  proposed  kinetic  theory 
based  model  for  in-tube  condensation.  The  second  objective  is  to  develop  an  energy-based  model 
for  in-tube  condensation.  The  focus  of  the  energy-based  model  is  the  determination  of  an 
expression  for  the  mass  velocity  Jm. 

Derivation  Of  An  Alternative  Heat  Transfer  Equation  For  Pure  Fluids  And  Near 
Azeotropes 

In  deriving  a practical  and  usable  expression  for  Jm,  the  Kraay  and  Rite  (1997)  theory  contains 
several  unrealistic  assumptions.  Consequently,  this  paper  proposes  an  alternative  method  for 
calculating  Jm.  The  new  method  retains  the  basic  approach  that  has  been  outlined  by  Kraay  and 
Rite’s  (1997):  the  total  heat  flux  is  a sum  of  the  single  phase  subcooling  and  the  two-phase 
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condensation.  Both  the  sensible  and  latent  heat  fluxes  are  based  on  an  enthalpy  difference,  rather 
than  a temperature  as  the  driving  potential  for  heat  transfer.  The  method  proposed  in  this  study 
differs  from  Kraay  and  Rite’s  (1997)  because  the  coefficient  7w’s  expression  is  derived  from 
mass  and  heat  balances  rather  than  from  the  kinetic  theory. 

As  stated  above,  the  total  heat  that  a condensing  fluid  exchanges  with  the  cooling  fluid  can  be 
separated  into  the  latent  and  the  sensible  heats: 

<l"=Jm„ifg  (3) 

and 

Jmsens  (if  -ijw)  (4) 

The  potential  for  the  latent  heat  transfer  is  the  enthalpy  difference  between  the  vapor  and  the 
liquid,  while  the  potential  for  the  subcooling  energy  transfer  is  the  enthalpy  difference  between 
saturated  liquid  and  subcooled  liquid  (wall  conditions).  The  corresponding  Jm ^ ^and  Jm  are 

functions  of  the  distance  along  the  pipe  (z),  because  these  coefficients  vary  with  the  vapor 
quality. 

The  expression  for  Jm { can  be  found  by  a simple  mass  balance.  Figure  1 represents  an  element 

of  liquid  of  length  dz,  which  lies  along  the  tube  wall.  Liquid  enters  the  element  from  the  left  side 
and  flows  along  the  wall.  Condensed  vapor  ( rhL[ ) enters  the  top  to  the  element  and  is  equivalent 
to  Jmi  ? multiplied  by  the  area  of  the  liquid-vapor  interface  ( dA ).  Liquid  exits  the  right  side  of  the 

element  at  a location  z+dz,  (i.e.,  rhL[+dz). 

The  mass  that  goes  from  vapor  to  liquid  per  unit  time  and  per  unit  area  (7  ),  can  then  be 

written  in  terms  of  the  variation  of  the  liquid  mass  flux  mL[  along  the  pipe  as: 

7 I -dA\  = dm,  I (5) 

Regardless  of  the  flow  pattern  of  the  fluid  inside  the  pipe,  the  area  of  the  interface  between  liquid 
and  vapor  can  be  written  as  2nr*dz. 

Jm  ■ 2 ■ 7t  ■ r*  -dz  = dm,  (6) 

mlut  L V ' 


7 


(8) 


The  relationship  between  the  vapor  quality  x and  the  liquid  mass  flux  is: 

dx  1 dmL 
dz  mwt  dz 

Substitution  of  Eq.  6 into  Eq.  8 gives  an  expression  for  J : 
rhtnt  dx 

J tot_ 

2 • K ■ r*  dz 


(9) 


The  relationship  between  the  latent  heat  transfer  coefficient  and  the  sensible  heat  transfer 
coefficient  is  found  by  means  of  an  energy  balance.  Consider  the  control  volume  shown  in  Fig.  2. 
The  enthalpy  of  the  bulk  fluid  is  assumed  to  be  a linear  combination  of  the  wall  and  saturated 
liquid  enthalpy,  as  shown  in  Eq.  10: 

i\z=ij\z  k+ifvlzfl-k)  (10) 

The  coefficient  k can  have  any  value  between  0 and  1,  while  ifw(z)  is  evaluated  at  the  wall 
temperature  and  at  a pressure  which  is  90  % of  the  pressure  of  the  vapor  phase.  This  last 
assumption  is  validated  by  the  fact  that  the  ratio  of  the  pressure  of  the  liquid  just  under  the  vapor- 
liquid  interface  to  that  of  the  gas  is  about  0.9  for  a low  velocity  of  the  vapor  toward  the  interface 
(Ytrehus  and  0stmo,  1996)  as  predicted  by  the  kinetic  theory  of  condensation. 


The  heat  balance  from  Fig.  2 has  the  expression: 

d(iiat  +mL[-i\+dmL[-if[  = mL[+dz-i[+dz+dqtot  (1 1) 


Because  dqtot=dqiat+dqsens,  Eq.  1 1 can  be  simplified  as: 
d[rtir\-i\  ] 

dqsen  s = Jz dz  + dmL[ -i\  ( 1 2) 

Dividing  by  dz  and  expanding  the  derivative  of  the  first  term  on  the  right  side,  Eq.  12  becomes: 


dq  sen.; 

dz 


d[<  |J. 

dz 


d[H  J 

dz 


mL  lz  + 


dmL  \z 
dz 


L I 


/ ‘z 


(13) 


Dividing  each  term  of  Eq.  13  by  the  circumference  of  the  vapor-liquid  interface,  the  left  side  of 
Eq.  13  represents  the  sensible  heat  flux  through  the  liquid-vapor  interface  area. 


(14) 


dqsens  1 _ l d[mL[]  1 d[i{ ] 

dz  2jrr*  27Tr*  dz  fz  z 2jtr*  dz  L " 
Since  Jm  is  defined  as: 

1 sen  s 


J = 


q sen  s 


•‘sens  / I _j  I 

lf'z  Lfw'z 


(15) 


Eq.  14  is  equivalent  to  the  following  expression  for  the  sensible  heat  molecular  flux: 


J = 


1 d[mL\ J 


27tr : 


dz 


iA~iK 
AA'z  ) 


mL\z  d[i\z] 


1 


27tr*  dz 


(16) 


To  correctly  sum  the  sensible  and  latent  heat  ( q”2nrdz ),  they  must  refer  to  the  same  area,  i.e.,  the 
area  of  the  liquid-vapor  interface  ( 2itr*dz ). 


q 


2tct  ■ dz 
2ttt  * -dz 


= q 


The  heat  flux  balance  is  then: 


1 


q" — = J if+ — 

■*  m,° ' fs  2nr*  dz  f 


1 4‘U, 


mL  K 


r mr*  dz  ' 2 7ir*  dz 

Substituting  Eq.  9 into  Eq.  18,  and  remembering  that  rhL\z  = mtnI  ■ (1  - x\z  ) gives: 


mtot  dx 


r * 2tct  * dz 
Simplifying  Eq.  19  gives: 


('A-'/O 


4rh.J * 


2tzt  * dz  dz 


m,. 


dM 


2ttt  * dz 


(1-x) 


2 Kr  dz 


dz 


(17) 


(18) 


(19) 


(20) 


Equation  20  represents  the  heat  balance  equation  for  the  condensing  fluid. 


The  data  presented  by  Kedzierski  and  Gonclaves  ( 1 999)  were  used  to  validate  the  model  given  by 
Eq.  20.  These  data  are  from  counterflow  tests  of  R32  flowing  in  a micro-fin  tube.  The  local  heat 
transfer  coefficient  was  measured  at  twelve  x-positions  along  the  tube  for  each  run.  Figure  3 
compares  the  experimental  heat  flux  versus  the  length  of  the  pipe  to  that  calculated  from  Eq.  20. 
The  value  of  k chosen  for  this  comparison  is  0.9.  The  agreement  between  the  measured  heat  flux 
and  the  model  chosen  for  the  heat  balance  is  within  2 %.  The  calculations  show  that  the  choice 
for  the  value  of  k in  Eq.  20  does  not  greatly  influence  the  results  ( q ” is  an  implicit  function  of  k, 
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because  the  bulk  enthalpy  i\z  is  a function  of  k).  This  is  due  to  the  fact  that  the  percentage  of  the 
sensible  heat  transfer  exchanged  over  the  total  (from  2 % to  7 %)  is  small  compared  to  the  latent. 
Therefore,  in  the  following  derivation,  the  value  of  k will  be  set  to  0.9. 


The  same  model  has  been  verified  for  other  refrigerants,  including  R134a,  and  for  a parallel  flow 
configuration  as  well.  For  all  the  tested  refrigerants  and  configurations,  the  agreement  between 
the  model  and  the  experimental  data  was  within  2 %. 


Equation  20  has  also  been  validated  for  zeotropic  refrigerant  mixtures  with  small  temperature 
glides  i.e.,  near-azeotropic  mixtures  such  as  R410A.  The  sensible  heat  due  to  the  change  in 
composition  of  the  vapor  has  to  be  added  for  mixtures  with  significant  temperature  glides,  e.g., 
R407C.  For  this  reason,  the  derivation  of  the  heat  balance  equation  for  the  zeotropic  refrigerant 
mixtures  is  dealt  with  in  the  following  section. 


Derivation  Of  An  Alternative  Heat  Transfer  Equation  For  Zeotropic  Mixtures 


During  condensation,  the  total  heat  released  by  a zeotropic  mixture  to  the  cooling  fluid  is  the  sum 
of:  (1)  the  latent  heat,  (2)  the  sensible  heat  in  the  liquid,  and  (3)  the  sensible  heat  due  to  the 
temperature  glide.  Once  the  variation  in  composition  as  a function  of  the  vapor  quality  along  the 
pipe  is  known,  the  enthalpy  of  the  liquid  can  be  written  as  function  of  the  vapor  quality  and 
inserted  in  the  heat  flux  balances  equation. 

Equation  14  can  be  rewritten  for  zeotropic  mixtures  as: 

1 


' K \if (x i, )-nx U) __J_ 


2 nr*  dz  " 27tr*  dz  u'  ' dz  2nr* 

After  applying  the  chain  rule  to  the  derivative  on  the  second  term  on  the  left  side  of  the  Eq.  21, 
the  final  expression  for  the  sensible  heat  transferred  in  the  liquid,  due  to  either  subcooling  or  to 
the  composition  change  (i.e.,  due  to  the  glide  in  the  liquid)  is: 


(21) 


1 dKU..  ........  1 d[i(x\)]dx  dq  1 

»-^r*  jz  —.mA  = 


27zr*  dz  1 ' z 1 27tr*  dx  dz  Lz  dz  27tr* 

To  find  an  expression  for  the  heat  due  to  the  glide,  a vapor  element  of  length  dz  is  to  be 
considered,  with  the  interface  included  in  the  element. 


(22) 
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For  a fixed  overall  mixture  composition,  performing  a heat  balance  on  the  vapor  element  gives: 


[ ;*'«  (A  )=  ^ Lft  dg  ( XlZ+dz  ) + d(l,at  + dmJ  f(A  ) + d(l  glide 

where  the  sensible  heat  due  to  the  glide  is  expressed  as  dqgude. 

Substituting  the  expression  for  dqiat  (Eq.  9)  into  Eq.  23  and  rearranging  gives: 

dd  glide 


= —rh. 


dis(  x\z  )dx 


dz  v dx  dz 

Since  rhv  I = rntot  ■ x\z , the  derivative  of  the  glide  heat  along  the  pipe  is: 


d<Z  cli, 


glide 


= ~mtoA 


dij  x\  )dx 


dz  z dx  dz 

For  zeotropic  mixtures,  the  heat  balance  equation  is: 

dq  tot=  dq  iat  + dq  sens  + dq  glide 
Substituting  into  Eq.  26  the  expressions  for  dq  ”iat  and  dq  ’ 
r 


and  the  above  term  gives: 


’<  dx 

2nr  * dz 


(ig(x l )~if(x[  ))■ 


K,  dx 

2nr  * dz 


(L(x[  )-i(x[  )- 


(23) 


(24) 


(25) 


(26) 


m di(  x[  ) dx  m dijxl  ) dx 

2nr*(  ^ dx  dz  27tr*  dx  dz 


(27) 


Combining  the  common  terms,  an  implicit  differential  equation  for  x that  applies  to  zeotropic 
mixtures  is: 


mu„  dx 

2k  r dz 


(i(x\z)-ig(x\z))  + x^ 


di(x  lz ) dig  (x  \z ) 


dx 


dx 


JJ 


= q 


(28) 


Description  Of  The  Experimental  Apparatus  Used  To  Validate  The  Enthalpy-Based  Heat 
Transfer  Model 

In  this  section,  the  enthalpy  based  energy  balances  for  pure  fluids  and  near  azeotropes  provides 
the  basis  for  correlation  of  the  local  tube-side  condensation  data  of  Kedzierski  and  Gonclaves 
(1999).  These  data  were  obtained  with  an  experimental  apparatus  used  to  establish  and  measure 
the  convective  condensation  for  R134a  and  the  near-azeotropic  refrigerant  mixture  R410A,  along 
with  its  pure  components  R32  and  R125. 
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The  Kedzierski  and  Gonclaves  (1999)  test  section  consisted  of  a pair  of  3.34  m long,  horizontal 
tubes  connected  by  a U-bend.  A magnetically  coupled  gear  pump  delivered  the  test  refrigerant  to 
the  entrance  of  the  test  section  with  a few  degrees  of  vapor  superheat.  Another  magnetically 
coupled  gear  pump  supplied  a steady  flow  of  water  to  the  test  section.  The  inlet  temperature  of 
the  water  loop  was  held  constant  for  each  test  with  a water-chilled  heat  exchanger  and  variable 
electric  heaters.  The  refrigerant  and  the  water  flow  rates  were  controlled  with  continuously 
variable  speed  pumps.  Redundant  flow  rate  measurements  were  made  with  Coriolis  and  turbine 
flow  meters  on  both  the  refrigerant  and  the  water-  sides. 

The  test  refrigerant  flowed  inside  a micro-fin  tube,  while  water  flowed  either  in  counterflow  or  in 
parallel  flow  to  the  refrigerant  in  the  annulus  that  surrounds  the  micro-fin  tube.  The  annulus  gap 
was  2.2  mm,  and  the  micro-fin  tube  wall  thickness  was  0.3  mm.  The  micro-fin  tube  had  60,  0.2 
mm  high  fins.  For  such  a geometry,  the  cross  sectional  area  is  60.8  mm2  giving  an  equivalent 
smooth  diameter  of  8.8  mm  and  a wetted  perimeter  of  44.23  mm.  The  root  diameter  of  the  micro- 
fin  tube  was  8.9  mm.  The  inside  surface  area  per  unit  length  of  the  tube  was  estimated  to  be  5.23 
mm.  The  ratio  of  the  wetted  perimenter  of  the  micro-fin  tube  to  that  of  a smooth  tube  with  an 
equivalent  hydraulic  diameter  was  1 .6.  The  fins  rifled  down  the  axis  of  the  tube  at  a helix  angle 
of  18°  with  respect  to  the  tube  axis. 

Heat  Transfer  Correlation 

Over  1 000  data  points  from  the  above-described  experimental  study  were  used  to  develop  a heat 
transfer  correlation  that  is  based  on  the  model  expressed  by  Eq.  20.  The  variable  to  be  correlated 

dx 

is  the  derivative  of  the  vapor  quality  with  respect  to  the  length  of  the  pipe,  — . The  quality 

dz 

. . dx 

derivative  represents  the  mass  of  the  vapor  condensed  per  unit  of  length.  The  values  of  — to  be 

dz 

correlated,  according  to  Eq.  35,  were  calculated  from  the  experimental  data  of  Kedzierski  and 
Gonclaves  (1999).  All  the  heat  exchanged  has  been  considered  to  be  latent.  Such  a choice  is 
mainly  based  on  the  observation  that  the  sensible  term  is  a small  percentage  of  the  total  heat 
exchanged.  From  the  analysis  of  the  data,  it  was  noted  that  for  the  pure  fluids  R32,  R134a,  and 
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R125,  the  sensible  heat  flux  corresponds  in  average  to  the  3 % of  the  total  heat  exchanged  in  the 
experimental  quality  range  (0.05-0.99).  For  mixtures,  the  fraction  of  sensible  heat  increases  to  6 

dx 

%.  Even  if  neglecting  the  sensible  heat  introduces  a systematic  error  in  estimating—,  the 

dz 

derivation  of  such  a quantity  from  the  actual  derivative  of  the  vapor  quality  introduces 
considerable  and  unpredictable  uncertainty  for  very  low  and  very  high  vapor  quality  values,  i.e.  at 
the  extremes  of  the  function  to  be  derived. 


The  actual  quantity  which  is  correlated  is  — ■ JaDh-  The  presence  of  Z)h  in  the  quantity  to  be 

dz 

dx 

correlated  allows  — • ja  ■ Dh  to  be  nondimensional.  The  Jakob  number  is  present  in  the  quantity 

dz 

dx 

to  be  correlated,  through  preliminary  analysis  of  D revealed  a variation  it  was  found  that  the 

dz  " 

exponent  on  Ja  that  was  nearly  Ja'1 . 


The  first  pass  for  the  search  of  the  form  of  the  correlation  is  done  using  only  those  data  points 
with  an  estimated  expanded  uncertainty  for  the  traditional  Nu  that  are  less  than  25%.  Next,  the 
686  omitted  measurements  are  predicted  with  the  same  regression  using  dimensionless 
parameters.  Those  excluded  measurements  that  are  predicted  to  within  ±25%  of  the  first 
correlation  are  reintroduced  into  the  regression  data  set.  New  coefficients  are  then  generated 
using  1317  measurements  of  the  regression  data  set.  The  hat  matrix  and  the  Cook’s  distance  are 
generated  for  the  second  regression  and  used  to  determine  the  data  that  has  have  high-leverage 
and  influence  on  the  model,  respectively.  The  diagonal  elements  of  the  hat  matrix  that  are 

greater  than  approximately  — , where  d is  the  number  of  model  parameters  and  n is  the  number 

n 

of  data  points,  are  designated  as  high  leverage  points  (Belsey  et  al.,  1982).  The  elements  of  the 
Cook’s  array  that  are  greater  than  — are  considered  to  be  influential  observations.  A total  of 

n 

1218  points  are  left  for  the  final  form  of  the  correlation. 
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The  form  of  the  correlation  was  obtained  by  following  the  approach  of  Cooper  (1984).  Cooper 
(1984)  suggests  that  the  fluid  properties  that  govern  nucleate  pool  boiling  can  be  well  represented 
by  a product  of  the  acentric  factor  and  other  dimensionless  variable  to  various  powers.  The 
present  correlation  has  the  form: 


Dh  - Ja  = a,  ■ Re*2  • Pr*3  ■ Sv*< 

dz 


f 

-Log  w 

\ 


W*5 


\PcJJ 


(29) 


where: 

a,  = 5.1081 
a2  = -0783242 
a3  = -0708893 
a4  = -0.606725 
a5  = 0.461457 

The  expanded  uncertainty  of  the  fit  for  95  % confidence  is  approximately  ± 21  %. 

dx 

Figure  5 compares  the  experimental  to  the  predicted  values  (Eq.  29).  The  expanded 


dx 

uncertainty  of  for  95  % confidence  was  approximately  ± 8 %.  For  high  values  of  the  quantity 

dz 


dx 

the  correlation  predicts  a different  behavior  for  counterflow  and  parallel  flow  data,  which  was 

dz 

also  the  case  in  the  work  presented  by  Kedzierski  and  Gonclaves  (1999).  The  difference  is 
appreciable  in  the  low  vapor  quality  regions  for  counterflow  conditions  and  in  the  high  quality 
region  for  parallel  flow  conditions,  as  shown  in  Fig.  6,  which  represents  the  trend  of  the 
derivative  of  the  vapor  quality  versus  the  vapor  quality  itself. 


Figure  7 compares  the  measured  condensation  Nusselt  number  with  the  Nusselt  number 
predicted  by  Eq.  29.  The  expanded  uncertainty  of  the  measured  Nu  for  95%  confidence  was 
approximately  16%.  Only  random  trends  are  observed  in  the  residuals  plots  against  each  of  the 
parameters  of  Eq.  29. 

The  present  correlation  was  compared  to  the  data  that  were  used  to  generate  it.  Figure  8 
compares  the  heat  transfer  coefficient  versus  quality  for  the  four  test  fluids  at  Gr  = 250  kg/m2s 
and  Tr=  40  °C.  The  lines  are  predictions  for  the  present  micro-fin  tube  geometry  which  were 
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obtained  form  the  correlation  of  the  data  given  as  Eq.  29.  In  general,  the  Nusselt  number 
decreases  for  decreasing  qualities  because  the  liquid  layer  becomes  thicker  as  the  quality 
decreases.  A greater  heat  transfer  resistance  occurs  for  the  thicker  condensate  films. 
Furthermore,  R32  exhibits  the  higher  heat  transfer  coefficient,  due  to  its  high  thermal 
conductivity  and  low  liquid  viscosity.  The  R125  has  the  poorest  heat  transfer,  while  the  mixture 
R410A  ( R32/R125,  50/50  composition  by  mass)  has  a performance  which  lies  in  between  those 
of  its  pure  components.  These  trends  are  consistent  with  those  reported  by  Kedzierski  and 
Gonclaves  (1999). 


From  their  relationship  to  the  traditional  heat  transfer  coefficient,  the  relation  between  the 


dx 

Nusselt  number  and  — is  expressed  by  Eq.  30: 

dz 

dx  pw-k,  1 

— = Nu 
dz 


m'cP  DhJa 


(30) 


Substituting  Eq.  30  into  Eq.  29,  an  expression  for  the  Nusselt  number  can  be  derived  as: 


D 

Nu  = — — • — • a,  ■ Re“2  ■ Pr*3  • 

Dh  Pwkt  ' 1 


-Log 


10 


( YY*5 
Pj_ 

<Pc  J, 


(31) 


mc„ 


The  hydraulic  diameter  cancels  out,  and  the  term  — — be  related  to  the  Peclet  number  (see 

Pj, 


Kedzierski  and  Gonclaves  (1999)): 
me 


P Nf-S  Nf-S 

P -Re-  Pr  ~ — =Pe-  * 


Pwkt  4cosa  • pw  ~~  4 cos  a ■ pw  4 

Substituting  Eq.  32  into  Eq.  30  the  final  expression  for  the  Nu  number  is: 


Nu  = — cc,  ■ Re“*  ■ Pra 3 • Svc 
4 


-Log 


10 


( Yf 

y Pc  )j 


where  the  new  exponents  for  the  Reynolds  and  the  Prandtl  numbers  are: 

a2  = -0.783242  + 1 = 0.216758 
0C3  = -0.708893  + 1 = 0.291107 


(32) 


(33) 
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Figure  9 represents  a comparison  of  the  correlation  expressed  by  Eq.  33  with  data  from  Eckels 
and  Tesene  (1998).  The  systematic  underprediction  of  the  Eckels  and  Tesene  (1998)  data  by  Eq. 
33  could  be  due  to  experimental  differences  between  the  data  that  was  used  to  generate  Eq.  33 
and  the  Eckels  and  Tesene  (1998)  data.  Baker’s  (1954)  flow  map  and  Soliman’s  (1982,  1983  & 
1986)  flow  transition  correlations  for  smooth  tubes  are  used  to  approximately  determine  the  flow 
conditions.  Baker’s  results  showed  that  79  %,  18  %,  2 %,  1 % of  the  data  are  respectively  in 
annular,  slug,  bubbly  and  wavy  flow  conditions.  Soliman’s  correlations  indicated  that  47  %,  30 
%,  and  23  % of  the  data  are  in  wavy  (stratified,  slug,  wavy),  annular,  and  mist  flow  conditions, 
respectively.  No  reduction  in  the  residual  standard  deviation  of  the  Nu  regression  is  obtained 
when  the  data  set  is  reduced  to  only  the  annular  flow  (as  indicated  by  the  Soliman  correlations) 
data. 

The  correlation  which  has  been  presented  in  Eq.  33  is  now  compared  with  the  one  proposed  by 
Cavallini  et  al.  (1997).  These  authors  correlate  the  Nusselt  number  with  5 variables,  Ree q,  Pr\,  Rx , 
Bo,  Fr.  These  are,  respectively,  the  equivalent  Reynolds  number,  the  liquid  Prandtl  number,  a 
geometrical  parameter,  the  Bond  Number  and  the  Froude  number.  For  an  accurate  description  of 
the  parameters  used,  refer  to  Cavallini  et  al.  (1997).  The  correlation  for  the  Nusselt  number  is: 

Nu  = 0.05  • Re£8  Pr1/3  Rx\Bo  ■ Fr )^26  (34) 

Figure  10  shows  that  if  the  above  expression  for  the  Nu  number,  given  by  Cavallini  et  al.  (1997), 

dx 

is  used  in  Eq.  30,  the  — values  are  overpredicted  by  approximately  200  % to  400  %. 

dz 
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Conclusions 

An  alternative  correlation  method  for  the  in-tube  condensation  heat  transfer  coefficient  which 
offers  some  advantages  over  the  traditional  method  has  been  presented.  The  alternative  approach 

dx 

correlates  a physical  quantity,  , rather  a mathematical  definition  as  is  done  in  the  traditional 

dz 

Newton's  law  of  cooling  approach.  The  correlation  that  results  from  the  alternative  method  is 
much  simpler  than  one  that  is  obtained  from  correlating  the  Nusselt  number  directly.  For 
example,  the  vapor  quality  does  not  appear  in  the  exponents  of  the  non-dimensional  parameters 
selected  for  the  correlation,  while  in  the  traditional  method,  the  presence  of  the  vapor  quality  in 
the  exponents  is  essential.  Another  improvement  with  respect  to  the  traditional  approach  is  that 
the  Jakob  number  does  not  appear  in  the  correlation.  The  benefit  that  is  realized  from  this  is  that 
an  iterative  procedure  is  not  necessary  to  calculate  the  temperature  difference  between  the 
saturated  vapor  and  the  tube  wall. 
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Figure  1-Mass  balance  on  a differential  liquid  element 


Figure  2-  Fleat  balance  on  a liquid  element 
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Figure  3-  Heat  flux  based  on  a heat  transfer  fluid  calometric  and  Eq.  (20) 
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dx/dz  [m'1]  predicted 


+dz  ig(^z+dz) 

> 


Figure  4-  Heat  balance  on  a vapor  element 


dx/dz  [m'1]  experimental 


Figure  5-  Comparison  of  measured  and  predicted  — , Micro-fin  tube,  R32,  R125,  R134a,  R410A 

dz 

condensation,  a=18°,  Dh=5.45,  Re=7715 
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dx/dz  [rrT1] 


Figure  6 
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Predicted  Nu 


Experimental  Nu 


Figure  7-Comparison  of  measured  and  predicted  convective  condensation  Nusselt  numbers. 
Micro-fin  tube,  R32,  R125,  R134a,  R410A  condensation,  a=18°,  Dh=5.45,  Re=7715 
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Figure  8-  Comparison  of  the  heat-transfer  coefficient  href  for  different  fluids  versus 

thermodynamic  quality 
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Nu  predicted 


Nu  experimental 

Figure  9-  Comparison  experimental  Nusselt  numbers  and  the  predictions  of  Eq.  (33). 
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predictions  of  Eq.  (34)  (Cavallini,  1997). 


Appendix  A:  Review  Of  The  Kinetic  Theory  Of  Condensation 


The  kinetic  theory  of  condensation  is  reported  in  several  publications  (Collier  and  Thome,  1 992, 
and  Wilhelm,  1962).  To  be  complete,  a brief  review  of  the  principal  concepts  of  the  kinetic 
theory  of  condensation  is  reported  below. 


Consider  a pure  saturated  vapor  with  pressure  Pg  and  temperature  Tg  in  equilibrium  with  its  own 
liquid  phase.  The  liquid  temperature  at  the  liquid- vapor  interface  is  7f.  Due  to  thermodynamic 
equilibrium,  the  liquid  moles  that  cross  the  liquid-vapor  interface  into  the  vapor  space  are 
balanced  by  an  equal  number  of  moles  that  move  from  the  vapor  into  the  liquid.  Condensation 
occurs  if  the  flux  of  molecules  from  the  vapor  to  the  liquid  phase  (Jm+)  exceeds  the  flux  of  the 
molecules  from  the  liquid  to  the  vapor  phase  ( Jm _).  Figure  A.l  illustrates  the  condensation 
process.  The  net  or  the  condensation  molecular  flux  across  the  vapor-liquid  interface  is: 

= (A.l) 

Using  the  kinetic  theory  of  condensation.  Hertz  (1882)-Knudsen  (1915)  demonstrated  that  the 
molecular  flow  rate  of  a stagnant  perfect  gas  through  an  imaginary  plane  is  given  by: 


M 

2tcR  ) 


\1/2 


T 1/2 


(A. 2) 


where  |7ra|  is  the  flux  of  mass  in  either  direction  perpendicular  to  the  plane,  M is  the  molecular 

weight;  R is  the  universal  gas  constant;  and  P and  T are  the  absolute  pressure  and  the  absolute 
temperature,  respectively. 


Substitution  of  the  Hertz-Knudsen  equation  (Eq.  A.l)  into  the  mass  flux  balance  (Eq.  A.l), 
where  Jm+  is  the  flux  from  vapor  to  liquid  and  Jm.  is  the  flux  from  liquid  to  vapor,  gives: 
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(A.3) 


Schrage  (1953)  modified  the  Hertz-Knudsen  equation  to  account  for  the  effect  of  a nonzero 
vapor  velocity  normal  to  the  interface  on  the  Maxwellian  velocity  distribution  by  inserting  a 
coefficient  (F).  Additionally,  two  more  coefficients  , <7C  and  Ge,  were  added  to  account  for  the 
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fraction  of  molecules  that  strike  the  liquid-vapor  interface  and  condense,  or  evaporate, 
respectively.  Schrage’s  version  of  Eq.  5,  therefore  is: 


J = 


( M ] 

1/2  " 

K 2tzRj 

r*.  P' 


'-pit 


-<y„ 


T, 


1/2 


(A.4) 


Values  for  F are  calculated  using: 

T - exp{-  a2)  + aK1/2[l  + erf  (a)]  (A.5) 

where  the  constant  “a”  represents  the  ratio  of  the  overall  speed  of  the  vapor  to  a characteristic 
molecular  velocity.  The  expression  that  relates  the  constant  “a”  to  the  molecular  flux  Jm  is: 


a = 


Pg(2M / RTg ) 


1/2 


(A.6) 
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Figure  A.  1 -Condensation  process 
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Appendix  B:  Kraay  And  Rite  Derivation 


To  derive  an  expression  for  Jm  for  use  in  Eq.  2,  Kraay  and  Rite  (1997)  simplified  Eq.  A.4  from 
the  kinetic  theory  of  condensation  by  making  several  assumptions.  Both  oc  and  oe  were  assumed 
to  be  equal  to  one.  This  assumption  is  commonly  made  due  to  the  difficulty  of  measuring  the 
values  of  the  condensation  and  evaporation  coefficients.  For  example,  Butterworth  (1977) 
suggests  that  setting  oc  equal  to  0.8  is  a good  engineering  approximation  for  all  circumstances. 


To  further  simplify  the  expression  for  Jm , the  temperature  jump  at  the  interface  was  taken  to  be 
equal  to  zero  (i.e.,  Tf=Tg ).  With  these  assumptions,  Eq.  A.4  becomes: 
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2P,~P * 

where  pg  is  the  actual  vapor  density,  and  pig  is  vapor  density  calculated  by  means  of  the  ideal  gas 
law. 

Because  the  kinetic  theory  does  not  provide  knowledge  of  the  pressure  at  the  liquid-vapor 
interface  from  the  liquid  and  vapor  sides  (P/and  Pg  respectively),  Kraay  and  Rite  (1997)  assumed 
that  Pf  was  equal  to  the  saturation  pressure  evaluated  at  the  wall  temperature  and  Pg  was  equal  to 
the  saturated  vapor  free  stream  pressure. 


Kraay  and  Rite’s  third  assumption  leads  to  a simplified  expression  which  describes  a 
condensation  process  where  the  ratio  of  the  vapor  mass  to  total  mass  (vapor  quality)  is  constant. 
However,  as  previously  stated,  in-tube  condensation  is  a strong  function  of  the  vapor  quality.  To 
remedy  this  inconsistency,  Kraay  and  Rite  ( 1 997)  assumed  that  the  relation  between  the  vapor 
quality  and  the  net  molecular  flux  was  linear.  At  a vapor  quality  of  100  %,  the  net  molecular  flux 
is  defined  by  Eq.  A.4.  Conversely,  at  a quality  of  zero,  the  net  molecular  flux  is  zero  because 
there  is  no  vapor  to  be  condensed.  The  resulting  equation  that  satisfies  the  above  two  quality-flux 
boundary  conditions  is: 


2Pp 


^ m,net  ^ m ^ 

2Pr  - PiB 


(B.2) 
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Equation  B.2  represents  the  molecular  flux  of  condensing  vapor  perpendicular  to  the  liquid- vapor 
interface.  The  molecular  flux  carries  with  it  the  energy  of  the  vapor  state  that  will  be  released  to 
the  vapor  upon  its  condensation  to  liquid.  The  heat  flux  due  to  the  phase  change  component  of 
the  condensation  process  can  be  expressed  as  the  product  of  Jm,net  and  the  latent  heat  of 
condensation: 


d(hat  _ j 2Pg 

dA  m 2 pg  - pig 


(B.3) 


The  sensible  heat  component  represents  the  amount  of  subcooling  experienced  by  the  condensed 
liquid.  The  sensible  heat  flux  can  be  expressed  as  the  heat  transfer  coefficient  times  the 
temperature  difference  between  the  liquid  at  the  wall  and  the  liquid  at  the  vapor-liquid  interface. 
The  expression  of  the  total  heat  flux  is  then  the  sum  of  the  phase  change  and  sensible  heat  flux: 


2Pg 

2Pg  ~ P‘g 


1 


(B.4) 


Superposing  a liquid  subcooling  heat  flux  to  a latent  heat  flux  allows  the  total  heat  flux  to  reduce 
to  the  only  the  sensible  heat  flux  at  the  cessation  of  condensation. 

Kraay  and  Rite  further  simplified  Eq.  B.4  by  coupling  the  latent  and  single  phase  terms  with  the 
following  assumption  (for  which  they  offer  no  justification): 


j 2p»  _ 

” 2Pi  - Pig  cr< 


Substituting  Eq.  B.5  into  Eq.  B.4  gives: 


(B.5) 
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— = J 2ps 
dA  "2ps-p,s 

Using  Eq.  B.6  Kraay  and  Rite  calculated  Jm  from  the  experimental  data  available  in  the  literature. 

Figures  B.l  and  B.2  compare  the  variation  of  the  Jm  and  /zr<?/ (traditional  heat  transfer  coefficient) 
coefficients  for  condensing  R32  normalized  by  the  values  at  the  test  section  entrance  for  parallel 
and  counterflow,  respectively.  The  data  are  taken  from  Kedzierski  and  Gonclaves  (1999)  and 
have  an  expanded  uncertainty  of  approximately  ±5  % for  95%  confidence.  These  data  show  a 
significant  variation  of  the  coefficient  Jm  with  heat  flux  along  the  pipe.  The  scatter  of  Jm  along 
the  pipe  is  comparable  to  the  variation  of  href.  It  was  believed  that  one  benefit  of  the  enthalpy- 
difference  approach  was  that  the  Jm  was  constant  with  respect  to  quality.  Figures  B 1 and  B.2 
illustrate  that  this  is  not  strictly  true.  Figure  B.3  also  shows  that  there  is  considerable  scatter  in 
the  variation  of  Jm  versus  the  vapor  quality  for  the  Kedzierski  and  Gonclaves  (1999)  data. 
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Figure  B.l -Normalized  Kraay  and  Rite  Jm  and  Normalized  h vs  the  distance  along  the  pipe 

Micro-fin  tube,  R32,  condensation,  a=18°,  Dh=5.45,  Re=7715,  Pr=0.278,  Parallel  flow 

configuration 
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Figure  B.2-Normalized  Kraay  and  Rite  Jm  and  Normalized  h vs  the  distance  along  the  pipe 
Micro-fin  tube,  R32,  condensation,  ot=18°,  Dh=5.45,  Re=7715,  Pr=0.278, 

counterflow  configuration 
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Figure  B.3-  Variation  of  Kraay  and  Rite’s  Jm  with  vapor  quality 
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Appendix  C:  Discussion  Of  Kraay  And  Rite  Assumptions 


In  Appendix  B,  the  assumption  made  by  Kraay  and  Rite  are  show.  In  this  Appendix,  the 
assumptions  made  in  Kraay  and  Rite  (1997)  are  discussed.  The  assumption  used  to  derive  Eq. 
B.6,  where  Pf  is  calculated  as  the  saturation  pressure  of  the  liquid  at  wall  temperature  conditions 
is  not  consistent  with  the  existence  of  subcooled  liquid  at  the  wall.  This  assumption  contradicts 
the  superposition  principle  used  to  obtain  Eq.  B.4,  where  the  liquid  at  the  wall  is  taken  to  be 
subcooled. 

Also,  an  inconsistency  arises  in  Eq.  B.4  when  the  configuration  of  the  heat  transfer  fluid  is 
counterflow.  For  parallel  flow,  the  total  heat  flux  diminishes  from  the  inlet  to  the  outlet  of  the 
pipe  due  to  a smaller  temperature  difference  between  the  heat  sink  and  condensing  vapor. 
Equation  B.4  correctly  reflects  this  trend  where  the  total  heat  flux  (left  term  in  Eq.  B.4) 
diminishes  with  increasing  distance  along  the  pipe  (z).  The  latent  heat  transfer  (first  term  of  the 
right  side  of  Eq.  B.4)  decreases  with  z due  to  the  linear  dependence  on  x (in  fact,  x diminishes 
along  to  pipe),  while  the  sensible  part  of  the  heat  flux  is  negligible  compared  to  the  total  heat 
flux.  Conversely,  for  counterflow,  the  total  heat  flux  increases  with  increasing  z,  while  the  latent 
heat  diminishes  with  z-  Consequently,  because  the  latent  heat  flux  typically  represents  over  90  % 
of  the  total,  Eq.  B.4  does  not  correctly  predict  the  trend  of  heat  flux  with  respect  to  pipe  length. 

The  other  observation  seen  in  Eq.  B.4  is  that  the  total  heat  flux,  as  well  as  the  sensible  heat  flux, 
passes  through  the  inner  area  of  the  refrigerant  pipe,  while  the  latent  heat  flux  passes  through  the 
liquid-vapor  interface.  Equation  B.4  does  not  account  for  the  difference  in  reference  surface  area 
for  the  latent  and  single  phase  heat  fluxes.  Consequently,  this  causes  the  energy  balance  not  to  be 
invalid. 
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